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Here we study properties of a homogeneous dipolar Bose-Einstein condensate in a weak anisotropic 
random potential with Lorentzian correlation at zero temperature. To this end we solve perturba- 
tively the Gross-Pitaevskii equation to second order in the random potential strength and obtain 
analytic results for the disorder ensemble averages of the condensate and the superfiuid depletion, 
the equation of state, and the sound velocity. For a pure contact interaction and a vanishing corre- 
lation length, we reproduce the seminal results of Huang and Meng, which were originally derived 
within a Bogoliubov theory around a disorder-averaged background field. For dipolar interaction 
and isotropic Lorentzian-correlated disorder, we obtain results which are qualitatively similar to the 
case of an isotropic Gaussian-correlated disorder. In the case of an anisotropic disorder, the physical 
observables show characteristic anisotropics which arise from the formation of fragmented dipolar 
condensates in the local minima of the disorder potential. 

PACS numbers: 03.75.Kk, 05.40.-a, 67.85.Hj 



I. INTRODUCTION 

Since the realization of Bose-Einstein condensates (BECs) in 1995 [11, there was a significant interest about 
effects of the disordered potentials on the properties of ultracold quantum gases. The reason for this is not only 
because of the unavoidable irregularities in the trapping potential induced by wire imperfections 0] , but also due 
to the fact that disorder can be generated and controlled using laser speckles Q . It is well known that cold atoms 
are a promising tool for simulating other physical systems 0] , in the sense of Feynman's quantum simulator [8]] . This 
applies also to the phenomenon of Anderson localization, which was originally used to microscopically describe the 
absence of diffusion in terms of disorder [9,]. It has a clear BEC analogue 10], which has been directly observed f6l. lTll . 
Also, localization inside BECs due to disorder created by atomic impurities on a lattice was studied theoretically [IJI 
and observed experimentally p^ . 

For a theoretical analysis of global dirty boson properties, different methods have been used to describe various 
limits, ranging from the Bogoliubov theory (l^ - [2l| to the Parisi replica method [23 - [26j . It turns out that long-range 
correlations within both the condensate and the superfiuid remain, despite the presence of disorder. However, both 
quantities are depleted due to the localization of fragmented condensates in the local disorder potential minima. For a 
strong enough disorder in a homogeneous system, the depletion increases to such an extent that even a critical disorder 
strength exists, above which a Bose-glass phase appears, consisting onl y of localized mini-condensates 1271-1321 . Effects 
of disorder were also studied for harmonically trapped BECs [H, [ll HI ll and BECs in optical latticesTmlllJg^^ 
W\ . while the temperature behavior of dirty boson properties was examined in Refs. [ij, [13, [12, [H, [ 



Realization of atomic dipolar BE Cs |39|-|41| with long-range anisotropic interaction has generated large interest in 
the theory of dipolar quantum gases |42|-[50j . Increase in the strength of dipolar interaction is possible by substituting 
atoms with magnetic dipoles by heteronuclear molecules, which have a strong electric dipolar moment in rovibrational 



ground state [5l| , or by inducing radiative coupling by placing dipoles into a resonator S^] . Dipolar condensates were 
studied in the case of isotropic disorder [s^, HI] , which yields characteristic anisotropies for both the superfiuid density 
and the sound velocity at zero temperature due to the anisotropy of the dipolar interaction. Typical disorder potentials 
realized in experiments are anisotropic, and this has been studied, to the best of our knowledge, only numerically and 
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for contact interaction [5J, |55| , whereas a 3d isotropic speckle disorder has also recently been proposed [56| . In this 
paper we analytically study the impact of a weak anisotropic disorder on a polarized dipolar BEC at zero temperature. 

To this end we proceed as follows. In SeclUwe calculate the lowest-order corrections of BEC properties due to the 
presence of disorder within a mean-field approach. For the sake of generality we consider an arbitrary two-particle 
interaction and a general disorder correlation function. In Sec. IIIII we specialize the developed formalism to dipolar 
interaction and a Lorentzian-correlated disorder in Fourier space. This yields for both the superfluid density and the 
sound velocity characteristic anisotropics, which should be measurable in an experiment. In Sec. II VI we present our 
conclusions and outlook for further related research. 



II. MEAN-FIELD APPROACH FOR WEAK DISORDER 

Bogoliubov quasiparticles and disorder induced fluctuations decouple in the lowest order [l3 - l2l| suggesting that 
disorder corrections can be calculated at zero temperature by neglecting quantum fluctuations and using a mean-field 
macroscopic wave function ip{r) governed by the time-independent Gross-Pitaevskii (GP) equation (ssl ]: 

-|^V2 + U{v) - M + y d?r'V{v - r')V'*(r')V'(r')) V-W = 0. (1) 

Here, m stands for the particle mass, ^ denotes the chemical potential, V{r — r') represents an arbitrary two-particle 
potential, while C/(r) describes an external disorder potential. Denoting the disorder ensemble average as (•), a 
homogeneous disordered system can be described, without any loss of generality, by a vanishing mean value {U{r)) = 
and an arbitrary correlation function {U{r)U{r')) = R{r — r'). With this the GP Eq. ([Ij is a stochastic nonlinear 
partial differential equation, where the statistics of the condensate wave function ipir) is governed by the statistics of 
the disorder potential U{r) [23|. Since ■0(r) describes the macroscopic occupation of the ground state, we assume it, 
without loss of generality, to be real. In addition to the statistical properties of the random potential we will assume 
that the macroscopic value of some physical quantity Pmac, obtained by coarse-graining of a microscopic quantity 
P(r), over a large volume V, gives the same result as the disorder ensemble average, namely: 

A„,ac = ^ / d'rA{r) = {A) . (2) 
'v 



V 



Here the length of the coarse-graining I ~ V^^"^ is assumed to be larger than both the correlation length a of the 
disorder potential C/(r), and the healing length ^ = ti^ / {2mng) , which represents the characteristic distance at 
which the condensate wave-function responds to some perturbation in the external potential: 

I » (3) 

In the definition of the healing length n represents the density of the fluid and g ~ Ai:Ti}aslm denotes the strength 
of a short-range interaction part of the two-particle interaction potential V(r — r') = g6(v — r') -t- . . . , expressed in 
terms of the s-wave scattering length . 

We consider the case of a sufficiently small random potential U (r) ^ gn, when the perturbative decomposition of 
the wave function of the system is justified: 

V'(r) =Vo(r)+Vi(r)+V2(r)4^... , (4) 

where corresponds to the correction of the wave function of order n in the disorder. Solving the GP equation 

([1]) in the zeroth order of U (r) gives 

whereas the first order correction is straight-forwardly calculated and its Fourier transform reads 

We note that its disorder ensemble average vanishes. Therefore, we also have to determine the second-order result, 
which turns out to be 

/• d^fc' £/(k - kOV^i(k') + V^o[2F(kO F(k)]^i(k')V>i(k - kQ 

In the following we use this systematic perturbative approach and calculate the respective properties of the dirty 
BEC. 
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A. Order parameter and condensate depletion 



In analogy to quantum field theory, the one-particle density matrix is defined as {tp{r)ip{r')) 32] . The macroscopic 
fluid density is the diagonal part of the one-particle density matrix, n = (^^ (r)) , whereas the condensate density is 
usually defined as the off-diagonal long-range order (ODLRO) parameter [32| 

no = lim {xp{r)ip{r')) . (8) 

|r— r'|— foo 

Coarse-graining of the one-particle density matrix (■0(r)-(/'(r')) over the fixed volume V before taking the limit does 
not change the result 

no= lim / d^nd^ra (V'(r + ri)?A(r' +r2)) . (9) 

|r-r'|-s.oo V Jv0V 

The integration commutes with the disorder ensemble average, and using Eq. ([2]), we obtain 

no = ((V(r)) i^Pir'))) = (^(r))^ . (10) 

The last equality follows from the fact that the average of an already averaged expression can be omitted. Therefore, 
the depletion of the condensate due to disorder, n — hq = (V'^) ~ (V')^ i is simply the variance of the wavefunction. 
Physically, this condensate depletion is due to the formation of fragmented condensates in the respective local minima 
of the random potential. Defining a separate Bose-glass order parameter by considering the ODLRO parameter of 
the two-particle density matrix 32] 

ino+qf= lim (V-WV^r')') = n2 , (11) 

|r— r' |— voo 

shows that the density of the fragmented condensates q defined in Eq. (fTTj) coincides with the condensate depletion 
n — no. To this end the disorder ensemble average is obtained along the same lines as Eqs. (|8|)- (fT0l) . Thus, we 
conclude that the localization phenomenon for quenched disorder follows already from a mean-field description of 
the dirty boson problem. Therefore, our mean-field approach represents a simplified derivation of the disorder- 
induced condensate depletion in comparison with the Bogoliubov theory of Refs. ^14i-i21i] . Note that disorder effects 
on Bogoliubov quasiparticles have recently been analyzed in Refs. [l9l, |2]|. 
The perturbative expansion (j4|) now yields for the particle density 

n = {^{rf) - + (Mrf) + 2^o (M^)) + ■ • ■ , (12) 

and, correspondingly, for the condensate density 

no = (V(r))' = i^l + (V'2(r)) + . . . . (13) 

This yields the condensate depletion 

n-no = {Mrf) + ... . (14) 
Using Eq. ([6|) we arrive to the following expression: 

[ d^k R{k) 

n — no — n / . ,„ ^ -|- . . . . (15) 

J (27r)3[Z|li + 2ny(k)]' 

Note that this represents a result for the condensate depletion in second order of the disorder potential for an arbitrary 
two-particle interaction potential and an arbitrary disorder correlation function. Specializing to the delta-correlated 
disorder i?(k) = R, and contact interaction ^(k) — g, Eq. ([T5l) reduces to 



which is the seminal result originally obtained by Huang and Meng within the Bogoliubov theory of dirty bosons 
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B. Equation of state 

Solving the equation {i^^ifJ-b)) — n^fit) for the chemical potential /i^ yields its dependence on the average fluid 
density nt, = Hb{n). We have introduced the notation nt,, denoting the "bare" chemical potential, because it diverges 
for uncorrelated disorder, regardless of the density n, as can be seen from inserting expressions ©-([T]) into the 
second-order correction (fT2|) : 

f,,^nVik = 0)- f ^ ^^^^^ ^ +.... (17) 

This unphysical ultraviolet divergence can be removed by renormalising the chemical potential [13] ■ If the density of 
the system vanishes, i.e. if there are no particles in the system, the energy needed for a particle to be added also has 
to vanish fi(n = 0) = 0. Therefore, we define the renormalized chemical potential according to 

^l{n)^^lb{n)~^lb{0). (18) 

With this we obtain, in second order of the disorder strength, the renormalized chemical potential: 

^^nV{k^Q)+An -————^^ + . . . . 19 

J (2^)3 rj^^n^^ 2ny(k)] 

which does not contain an ultraviolet divergence. 

For calculating the sound velocity later on we will also need the expression for the compressibility of the fluid, or 
its inverse given by dfi/dn. Note that from Eq. (|18|) it follows that the obtained result does not depend on whether 
we use fj, or lib- Thus, from the perturbative expansion (I19|) we read off: 



^.y(k^o) + 4/^^&)^ 



y(k = 0)+4 ^J" .3 +■■■ ■ (20) 



C. Superfluidity 



Without disorder and at T = 0, the whole system is in a superfiuid state, moving with an arbitrary wavevector 
ks, which corresponds to the superfiuid velocity V5 = hks/m. By introducing disorder that moves with the velocity 
hku/m, some part of the fiuid will be moving together with it. The normal, i.e. non-superfluid, component of the 
fluid UN is defined as the part that moves together with the disorder, while the superfiuid component ns is defined 
as the fraction of the fluid that moves with the superfiuid wavevector kg. Therefore, the macroscopic current density 
(j(r)) can be separated in this two-fluid picture as follows: 



(j(i")) = nsks + UNki 



(21) 



The averaged current density (j(r)) can be obtained by analyzing the underlying time-dependent GP equation for the 
system: 



2m 



U r 



ku^tj + J d'r'V{r~r')%ir\t)^s{r',t) 



dt 



(22) 



where the condensate wave function VPg is a product of some as yet unknown function ipg and a plane wave with 
wavevector ks that corresponds to the clean-case solution: 



*5(r,i) =e'''^>5(r,t) 



(23) 



Substituting the ansatz 
to 



into Eq. (1^^ . changing variables via x = r — kij—t and introducing K. = ks — kjj leads 



- i—KV + C/(x) - Ai + / d'^x'Vix - x')Vs(x')V's(x') 
2m m J 



V's(x) =0. 



(24) 
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Although 4's should in general depend on t, it can be shown via mathematical induction on the perturbative solution 
that all orders of V's(x, turn out to be time- independent [s^l- Note that -05 does not depend explicitly on the 
wavevectors kg and kjj, but only on their difference K. Here, we are only interested in small values of K and, 
therefore, expand ips ips — ip + pK + . . . , with p — (9V's/9K)k=o- An explicit equation for p can be obtained by 
performing the derivative of Eq. (1241) with respect to K , yielding 



-|;^V2p(x) - i^V^(x) + [{/(x) - ^l] p(x) 

+ / d^x'V{x - x') { [p* (x') + p(x')] i^i^'M^) + V^(x')2p(x)} = . (25) 



transforms into 



If we take into account Eq. (j23l) . the standard definition of the current density 

(j> = ^ (*sV*s - ^sVn) (26) 

(j) - mi^s) ks + ^ (^sVV^s - i^s^rs) , (27) 

which then can be further reduced to 

(j) = nks + ((■0V (g) Imp) - (V0 <g) Imp))K + . . . . (28) 

In the last line we have neglected higher than linear orders in hjj and kg. 

For small disorder strengths, if we expand Eq. ((28)) with respect to U up to second order, take into account the 
homogeneity of our problem that leads to di (p2(x)) — and note that in zeroth order tps does not depend on K, 
thus leading to po — 0, we obtain: 

(j) = nks + {{tpiV ® Impi) - {VtPi ® Impi))K + . . . , (29) 



where also p is expanded in the disorder strength as p = po + pi + P2 + . • . • Solving the imaginary part of Eq. ([25)) 
in the first order in U gives the Fourier transform of the Impi: 

(Impi)(k) = 2i^0i(k). (30) 



Thus, together with the solution for -0i given by Eq. (|6]), we obtain from Eq. (|29)) and comparison with Eq. (|2ip the 
normal fiuid density in the form 

, f d^k k^k R(k) , , 

nN = ^n ——:r—- — ^ + .... (31) 

J (2^)3 fc2 ^2^^2ny(k)]' 

Note that in general the non-superfiuid component is represented by a tensor (ssj . 

In the case of a cylindrically symmetric system, we can choose the symmetry axis as the z-axis and denote the 
polar and the azimuth angle by 9 and if, so integrating Eq. ()3ip in spherical coordinates with respect to if yields the 
angle dependence 



sm 



„2ir /■2n / sin^ 9 cos^ if sin^ 9 sin tp cos </? sin 9 cos 9 cos ip 

9 dtp eke J = sm6 dip \ sin^ 9 sin p cos ip sin^ 9 sin^ p sin 9 cos 9 sin p 
•^0 \ sin 6* cos cos (/3 sin cos 6' sin (/s cos^ 6* 



/ 7r(l - cos2 e) \ 

= sin6' 7r(l-cos2 6i) . (32) 

\ 27rcos2 6' / 

If both V{k) and i?(k) are 6'-independent, i.e. if we have spherical symmetry, integrating Eq. ()32p with respect to 9 
leads to a solution in second order of the disorder potential [l3 - l2Tj | : 

4 

fi-N = -^{n-no)I ■ (33) 

This result shows that the superfluid depletion will be by a factor of 4/3 larger than the condensate depletion. Thus, 
the localized fragmented part of the fluid hinders the superfluid to move. 
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D. Sound velocity 

In the mean-field approach, we can also define the sound velocity by perturbing the time- independent solution with a 
small time-dependent variation. It is expected that sound waves with wavelengths of the order of the correlation length 
would scatter and interfere due to disorder hills and valleys, making the sound velocity impossible to define precisely. 
Locally, the sound waves would have the same speed as in the clean case. For sound waves with wavelengths much 
larger than the disorder correlation length, the sound velocity can be calculated using the hydrodynamical approach 
[5 Si] . Hydrodynamic equations are valid in the macroscopic regime and can be used only for slowly varying quantities 
that do not depend on the specific microscopic realization. Spatial averaging over distances much larger than the 
correlation length and much smaller than the wavelength solves the problem. Assuming that it gives the same result 
as the disorder ensemble average, we obtain the hydrodynamic equations for the macroscopic, i.e. disorder averaged, 
quantities in the form 

^^^ + V(ns(x,t)vs(x,t)) = 0, (34) 

m — -I- VI ^ +/i(n(x,t))l =0, (35) 

where n denotes the macroscopic density, and the disorder velocity k[/ is taken to be zero. If we write densities and 
the superfluid velocity as sums of homogeneous equilibrium values and small variations, 

n(x, i) = n + (5n(x, i), (36) 
ns(x,t) = ns + (5ns(x,t), (37) 
vs(x,t) = (5vs(x,t), (38) 

as well as neglect second-order terms in the variations, we get the following linearized system of equations: 

d6n{x, t) 



dt 



VihsSvs{^,t)) = 0, (39) 



(3Jvs'(x, i) 1 , „ , 1 du „ , , , , 

= Vfi{n + Sn{^,t)) = J-VSni^,t). 40 

at m mon 

Taking the time derivative of Eq. (I39p and substituting the expression for the superfluid velocity variation from 
Eq. (j40p we obtain the generalized wave equation 

^^^-^^V(n,VMx,t))^0. (41) 
m on 

From the above equation we deduce that the sound velocity in the direction of on unit vector q is given by 

cl = -^ci'nsci, (42) 
^ m on 

where the tensorial property of the superfluid density has been taken into account. In order to further evaluate 
the sound velocity (H^ for small disorder, the perturbative results for both the inverse compressibility (PH)) and the 
superfluid density following from ([3T|) have to be taken into account. 

III. DIPOLAR INTERACTION AND LORENTZ-CORRELATED DISORDER 

In this section we will specialize the previously developed perturbative formalism and consider BEC systems in the 
presence of two different anisotropics, namely an anisotropic dipolar interaction between the analyzed particles and an 
anisotropic disorder potential. The latter is widely studied and physically motivated, for instance, by the anisotropy 
of the laser-speckle potential [1, Q. In order to obtain analytical results, we model the disorder correlation function 
by a cylindrically-symmetric Lorentzian in a Fourier space 

m- . , J . (43) 



The lengths ap and az are respectively the perpendicular and the parallel correlation length. This function is not phys- 
ically realistic, but the corresponding results qualitatively coincide with the case of an isotropic Gaussian-correlated 
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disorder, which was numerically calculated in Ref. [53|, and we expect that all phenomena that appear here would 
also appear qualitatively for a true laser-speckle correlation function in a setup where it decays monotonously with 
distance. 

Assuming that the van der Waals forces between the atoms can be approximated at low energies by an effective 
contact interaction, the interaction potential in the presence of an external field that aligns the dipoles in a direction 
m takes the form [59| 

F(r) = g(5(r) + -^(l-3cos<^(m,r)) , (44) 

where 4>{m, r) represents the angle between vectors m and r , and C^d denotes the dipole-dipole interaction strength. 
In the case of magnetic dipoles Cdd = /^o"^^, with /ip being the magnetic permeability and the magnetic dipole 
moment m, whereas for electric dipoles we have Cdd = <P / (-a, with the vacuum permeability eg and the electric dipole 
moment d. Introducing the ratio of the dipole-dipole and the contact interaction e = Cdd/35 , and taking the Fourier 
transform of the potential, we obtain 

F(k) = g { 1 + e [3 cos^ k) - l] } . (45) 

The Huang and Meng result [3| for the condensate depletion (|T6l) is linear in R, and therefore we will compare the 
relative change of physical quantities due to disorder to the relative change of the condensate density. To this end we 
define a dimensionless disorder correction of the quantity A as 

— - Arf 

Aa = hm ^ . (46) 

n 

Here stands for the quantity A in the clean, isotropic system, while Ad denotes a possible dimensionless anisotropy 
factor due to the dipolar interaction and tihm represents the Huang-Meng depletion from Eq. (|16l) . Corrections defined 
in this way are expressed in terms of only three parameters: the relative dipole-dipole interaction strength e, and the 
correlation lengths in units of the healing lengths, Zp,z = y/2ap^z/S,- We consider systems with an overall cylindrical 
symmetry, where the disorder symmetry axis is parallel to the direction of dipoles. Otherwise the angle between them 
would be a fourth parameter that would have to be taken into account. 

The general case can be calculated analytically for all observables of interest, but the results are too cumbersome 
to be displayed here [131 . In the following we restrict the results and their discussion to two special cases, namely the 
pure contact interaction with anisotropic disorder, and the dipolar interaction with isotropic disorder. The general 
situation of having both dipolar interaction and anisotropic disorder just leads to quantitative but not qualitative 
changes. 



A. Condensate depletion 

We now use the disorder correlation function and the interaction potential defined by Eqs. (I43p and (j45p in order to 
calculate the disorder corrections of the condensate density, chemical potential, inverse compressibility, and superfluid 
density, given by Eqs. (|15l) . p9|) . (|20)) . and ((3T|) . Taking into account Eqs. (fT6|) and p6)) . making a substitution 
k — k^/\/2, denoting the direction of the disorder symmetry by d and introducing direction-dependent anisotropy 
functions r and v 



I zf, sin^ 0(d, k) + zl cos^ <;()(d, k) , (47) 
V = Vl + e(3cos2(/)(m, k) - 1) , (48) 
Eq. (fTSj) yields the dimensionless disorder correction for the condensate depletion in the second order 

. _ ( ^ 
-Ano-SnJ (2^)3 (k2+ ^2)2(1 + k2^2)- (49) 

Assuming that the direction of dipoles is parallel to the disorder symmetry, i.e. d||m, yields that the whole system is 
cylindrically symmetric. Writing Eq. (|49)) in spherical coordinates {k,6,ip), integrating with respect to k and (p and 
changing the variable t — cos 9, leads to: 

- A„, - / dt— ^ (50) 

Jo (1 + vr) 



8 





FIG. 1: (Color online) Condensate depletion due to weak disorder for: (a) anisotropic disorder and contact interaction from 
Eq. (|53p and (b) isotropic disorder and dipolar interaction from Eq. (|54p . 



with functions r and v from Eqs. (j47p and (|48p having the new form 



V = Vl~ 



3et2, 



(51) 
(52) 



The two special cases of contact interaction {v = 1) and isotropic disorder {r = z) can be solved by using Euler 
substitutions r = xt + Zp and v = xt + y/l — e, respectively, in Eq. (j50p . which leads to an integral of a rational 
function with respect to x. Analytic results in the two cases are given by: 



^"o(e-O) 



(Zp-l)(Zp + l) 



2z2 



(Zp + l){Zp + Zz) 



T 



1 z. 



'^riQ {Zp 



z{l-\) 



Zp -\- \ Zz ~\~ Zp 



1 



T 



Zz + 1 

z8-\ 



(-l + z2,52)(l- A + z5(l + A)) (5(-l + z<5)(l + z(5)2 + 1 



with (5 and A introduced by (^2 = 1 — e and e = 



4A 



l-2A+3A^ 

new function well defined for positive values and analytically continuable for 



to simplify the expressions, and T{x) = 
1< a; < 0. 



arc tan 



(53) 
(54) 
being a 



In Fig. [T] we have displayed the results for the disorder correction of the condensate depletion ([55)) and Note 
that for e = and z = the result is 1, meaning that it coincides with the Huang and Meng result. For increasing 
correlation length the depletion decreases, and vanishes for infinite correlation length, as expected, since then the 
disorder is flat. The tiny asymmetry in Fig. [T] a) comes from the fact that Zp describes two spatial dimensions 
and, therefore, has a more pronounced effect on the depletion than Zz which represents only one spatial dimension. 
Increasing the relative dipole-dipole interaction leads to a larger depletion and, eventually, when it reaches the same 
order as the contact interaction, the BEC collapses, which is seen in our model as a divergence in Fig. [T]b). Note 
that the condensate depletion can be measured in matter interference experiments, where the fragmented part of the 
fluid contributes with a random phase and, therefore, reduces the contrast of the interference pattern. 



B. Equation of state 

Corrections of the chemical potential and the inverse compressibility can be calculated from Eqs. and ([20]) in 
the same way as the condensate depletion. The analytic results for the correction of the chemical potential are given 

by 

A.(. ^ 0) ^ ,T (^^) + f , (55) 

[Z^ ~ l)(Zp + l)[Zp + Zz) \Zp + l Zz+ ZpJ Zp + Zz \Zp + ZzJ 

_ ^ -1 + A + 2zd{-l - A + zS[l ~ A + zSjl + A)]} 

Ap[zp -zz-z) - 2 ^2^2) [i^x + z6il + A)] 

'^''^'^n-X} + ^,,^tl±^T(^x] . (56) 

^ ' z^S{-l + zd){l + zS)^ V^-^+l 



The disorder corrections of the chemical potential are shown in Fig. [2j The correction increases with increasing 
disorder, regardless of the strength of the dipolar interaction and disorder correlation lengths. This is due to the 
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FIG. 2: (Color online) Change in the chemical potential due to weak disorder for: (a) contact interaction and anisotropic 
disorder from Eq. (|55|) and (b) isotropic disorder and dipolar interaction from Eq. (|56|) . 



repulsive interparticle interaction, which has a higher potential energy when the fluid is less uniform. The correction 
of the chemical potential in the case of the contact interaction, shown in Fig. [2] a) , has a similar dependence on the 
correlation lengths as the condensate depletion in Fig.[T]a), while, according to Fig.[5]b), the dipole-dipole interaction 
does not have a significant effect. This is due to the fact that the dipolar interaction contributes partially as attractive 
and partially as repulsive, thus leading only to a small net effect. Note that the chemical potential in the clean case 
is anisotropic as can be seen from Eq. (1191) and the directional dependence of the limit k rightarrowO in Eq. (|45l) . 
This is discussed in more detail in Ref. [50|. 

The corresponding results for the inverse compressibility are: 



Ae^(e = 0) 



A a,i (zp = Zz = z) = 

On 



y2\ 1 



3 + (2 + 



+ 



-4 + ^^) 



(-1 + A) {1 - A + z^S^\2 - 2A + Zz6{\ + A)]} 



-T 



1 Zr, 



Zp + 1 Zz + Zp 



z{-l + z^S^f[l-X + z6{l + X)]'^ 
They represent intermediate results for calculating the sound velocity in subsection IIII D[ according to Eq. (|42| 



(57) 
(58) 



C. Superfluidity 

The disorder correction of the superfluid density tensor can be separated into a perpendicular and a parallel 
component after integration with respect to ip, using Eq. (j32p we get 



A. 



















2A„o f 


1 


• 






-1 


•) 


V 

















(59) 



where A„(, is already calculated in Eqs. ([53| and (|54)) . and Isd is a new integral of the form 

t2 



dt- 

V {1 + vr)' 

which yields in the two studied special cases the following solutions: 



/sd(e = 0) 



1 



{-Zz + Zp) (Zz + Zp) 



-2 + zl 



Isd{Zp = Zz = z) = 



(1 + Zp) {Zz + Zp)' 

(-1 + A)3 



-T 



2 + Zz 

1 + Zz 



Zp -\- \ -\- Zp 



-T 



Zp -\- 



(-1 + A)= 
4z(52A[l - A + z(5(l + A)] ^z'^b^X 



T(-A) + 



(-1 + A)^ rj.(z5-^ 



-A 



(60) 



(61) 
(62) 



In an arbitrary direction of an unit vector q the superfluid density can be calculated using the tensorial superfluid 
density by ns(q) — (\hs(\. In the case of the cylindrical symmetry it is reduced to 



ns(q) = ns„ sin^ 0(q, e^) + n^, cos^ 0(q, Gz). 



(63) 
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FIG. 3: (Color online) Change in the (a), (c) perpendicular and (b), (d) parallel superfluid density due to weak disorder for (a), 
(b) contact interaction and anisotropic disorder and (c), (d) dipolar interaction and isotropic disorder. Red lines show values 
of parameters where the superfluid depletion becomes equal to the condensate depletion. 




FIG. 4: (Color online) Critical values of the parameter e for (a) perpendicular and (b) parallel superfluid density, at which the 
superfluid depletion becomes equal to the condensate depletion. 

Thus, obtaining the disorder corrections for ns^ and ns^ is sufficient for recovering thesuperfluid depletion in any 
direction. From Eq. (|59|) we directly read off 

-As, = -2A„„-2/,d, (64) 
-As, = 2/,d. (65) 

Note that the "-" sign in front of As, , and A„q suggest that the change of the superfluid densities and the condensate 
density are negative, or, equivalently that the depletion is positive. For isotropic systems Ag, and Ag, are both equal 
to I A„o, as can be seen from Eq. (1331) . Due to an anisotropy, however, there is a range of correlation lengths and 
relative dipolar interaction strength, where the superfluid depletion is smaller than the condensate depletion. Some 
particles from the fragmented fluid contribute to superfluidity, suggesting that they are not localized indefinitely, but 
have some finite localization time, which was calculated in Ref. [32J within the Hartree-Fock theory of dirty bosons. 
The superfluid depletion in the case of contact interaction and anisotropic disorder shows a similar behaviour as 
the condensate depletion, as can be seen in Figs. [3] a) and[3]b). For isotropic disorder and dipolar interaction the 
depletion of the perpendicular component is similar to the condensate depletion Fig. |3]c), but the depletion of the 
parallel component decreases as the relative interaction strength increases, as shown in Fig. |3]d). The red lines in 
Fig. [3] show where the superfluid depletion becomes equal to the condensate depletion. Fig. |3] shows the values of e 
for which this change happens. 

Although defined only for systems without a trap, the above calculated superfluid density could be extended to 
the trapped case in the simplest way by assuming that it depends only on the local density. If we turn on a slowly 
moving disorder for a short time t, such that vt is much smaller than the size of the trap, before switching off the 
trap, this would change the momentum distribution which could be afterwards reconstructed from a time-of-flight 
measurement. In this way our predictions for the superfluid density in such a system might become observable in 
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FIG. 5: (Color online) Sound velocity for dipolar interaction 
disorder. 




FIG. 6: (Color online) Corrections of (a) perpendicular and 
contact interaction. The red lines show values of the correlatioi 
correction vanishes. 




(a) clean case and (b) correction due to weak delta-correlated 




(b) parallel sound velocity for weak anisotropic disorder and 
lengths and the relative dipole-dipole interaction for which the 



experiment. 



D. Sound velocity 



The corresponding disorder correction of the sound velocity in the direction of the unit vector q can be calculated 
using Eqs. ^ and (|46l): 



gn 
m 



K,(0) nHM / ^q(O) T/C ^ , . 

1 q A„„q + Aa^ 

n \ g a 
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which leads to 



g 1^ 



(66) 



(67) 



where Vq(0) — limfc^o y{kc\) denotes the directional dependence of the potential V on q, according to Eq. (|45t . The 
anisotropy factor due to the dipolar interaction = Vf^{Q)/g and the delta-correlated disorder correction (53| are 
plotted mFig.m 

The anisotropy of disorder comes into play in a simple way. From Eqs. (j42p and (|63p it follows that the sound 
velocity can also be separated into a parallel and a perpendicular component 



with 



(q) = ci sin^ ^(q, ) + cos^ (/)(q, ) 



1 5/i 

m on 



and the corresponding dimensionless disorder correction 

A,2 =A8m 



(68) 



(69) 



(70) 



By introducing a weak disorder, the sound velocity changes via two competing effects: the decrease of the compress- 
ibility, i.e. the increase in the inverse compressibility, from Eq. (j20p . which tends to increase the sound velocity, and 



12 



the decrease of the superfluid density (negative value of Ag^ ^ ), which tends to decrease the sound velocity. The results 
are shown in Fig. [51 The red line shows where the second order sound velocity correction vanishes. For small disorder 
correlation lengths the decrease in compressibility is dominant. These corrections can be experimentally measured by 
determining the phonon dispersion relation by using Bragg spectroscopy 0, l60l , Isij . 

IV. CONCLUSIONS 

We have analyzed in detail how the anisotropy of both the dipolar interaction and the presence of disorder affects the 
directional dependence of different physical observables of dirty Bose-Einstein condensates. Using at zero temperature 
the mean-field approach, we have calculated the condensate depletion due to disorder, as well as the corresponding 
corrections to the equation of state, the superfluid density, and the sound velocity. In particular, we have discussed 
the consequences for the sup erfluid density, which becomes a tensorial quantity as a linear response to the moving 
disorder. Whereas Ref. [53[ analyzed a dipolar BEC in isotropic disorder potential, we have shown here that the 
anisotropic disorder provides a separate origin for the tensorial nature of the superfluid density. We have found that 
a large enough disorder anisotropy can even make both the parallel and perpendicular superfluid density component 
larger than the corresponding condensate density, which happens in the case of dipolar interaction and isotropic 
disorder only for the parallel component [s^. These initial results necessitate further studies, as they contribute 
to the overall physical picture in which the localization of bosons in the respective minima of the disorder potential 
occurs at a characteristic time scale (32j . This localization time remains to be analyzed in more detail in a forthcoming 
publication. 
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